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Abstract. We model the formation of high-mass stars, specifying the 
! accretion rate in terms of the instantaneous and final mass of the star, the 

ambient pressure of the star-forming region and the form of polytropic 
pressure support of the pre-stellar gas core. The high pressures typical of 
' Galactic regions of massive star formation allow a 100 Mq star to form 

, in ~ 10^ yr with a final accretion rate ~ 10"'^ Mq yr~^. By modeling 

' protostellar evolution we predict the properties of several nearby massive 

. protostars. We model cluster formation by applying this theory to many 

I stars. We use the observed intensity of outflows from protoclusters to 

. estimate the star formation rate, finding that clusters take at least several 

free-fall times to form; for a cluster similar to the Orion Nebula Cluster, 
Q^. we predict a formation timescale ~ 1 x 10^ yr. 

o 
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1. Introduction 

If there were no massive (m* > 8Mq) stars, the heavens would appear quite 
dull. Deprived of many of the heavy elements, it is not clear that Earth-like 
\ astronomers could even exist to contemplate the obscurity of such a universe. 

Fortunately, back in our own cosmos, we can be thankful for the existence of 
high-mass stars and perhaps devote some effort to understand them and how 
they come about. All massive stars appear to form in clusters, where they are 
greatly outnumbered by their lower-mass cousins, and so the questions of how an 
individual massive star forms and how the cluster forms are intimately entwined. 

In this article we first outline a theory for the formation of individual stars, 
including massive ones, by accretion (§2, taken from McKee Sz Tan 2002; here- 
after MT02). High pressures in the star- forming region result in rapid accretion 
rates, which allow massive stars to form in spite of the fierce glare of their radia- 
tion (Wolfire & Cassinelli 1987), and to do so swiftly before stellar evolution has 
time to age the newborn stars. In certain regions of the Galaxy, the required high 
pressures are known to result from the weight of hundreds to thousands of solar 
masses of strongly self-gravitating and turbulent gas, condensed into clumps of 
order a parsec across. A large fraction of this material inevitably collapses into a 
stellar cluster, with most of the mass partitioned into relatively low-mass stars. 
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The massive stars are expected to form preferentially at the center of the clump, 
where the pressures are highest. We apply our theory for individual star forma- 
tion to model the birth of the cluster (§§3 & 4, taken from Tan & McKee 2002a; 
hereafter TM02). Extrapolating models of outflows from low-mass protostars to 
those from higher masses, we use observations of protocluster outflow intensities 
to estimate the star formation rates and cluster formation times. 



2. Collapse of Turbulent Cloud Cores and Massive Star Formation 

Our basic assumption is that a star forms from the gravitational collapse of 
a molecular cloud core that begins in approximate hydrostatic equilibrium; in 
particular, it does not form by the coalescence of low-mass stars (Bonnell, Bate, 
&: Zinnecker 1998) or by the sweeping up of ambient clump material (Bonnell 
et al. 1997). On dimensional grounds we then expect that the protostellar 
accretion rate is given by 

= 0* — , (1) 

where tg = {2>tt /?>2GpY^'^ is the free-fall time and (/>* is a dimensionless con- 
stant of order unity. This equation has the same dependence on dimensional 
parameters as that for isothermal collapse, ~ ^j^h/^ (Shu 1977), if the 
thermal sound speed, Cth, is replaced by the virial velocity (Gm*/i?)^/^, since 
(Gm*/i?)3/2/G oc m=,(Gm*/i?3)V2 qc m*/tff. We have shown (MT02) that if 
the collapse is spherical and self-similar, then 0^, is quite close to unity provided 
that the value of p entering % is evaluated at the radius in the initial cloud that 
just encloses the gas that goes into the star when its mass is m*. We summarize 
these results below. Collapse to the star will naturally proceed via a disk, where 
we assume instabilities drive high inflow rates if the disk mass grows to be a 
significant fraction of the star's (Shu et al. 1990). 

First we allow for the disruptive effect of protostellar outflows on the gas 
core so that only a fraction ecore = m^/M of the core mass can accrete onto the 
star. The value of ecore is uncertain for massive stars, since it depends on the 
unknown properties of their protostellar winds. We assume that it is constant 
over the accretion history of the star, and for numerical estimates set ecore = 0.5, 
similar to the value calculated for low-mass stars (Matzner & McKee 2000). 

For the core, we assume a spherical power-law ambient medium, which 
ensures accretion is self-similar. With p oc r~^p and P oc r'^^p^ it follows that 
the core is a polytrope with P oc p'''p. In hydrostatic equilibrium (HSE) kp = 
2/(2 - 7p) and kp = -/pkp = 2jp/{2 - jp) (McLaughlin & Pudritz 1996). Let 
c = (Pj p)^!"^ be the effective sound speed. The equation of HSE then gives 
M = {kpcPr)/G and p = AcP /{2TTGr'^), with A = 7p(4-37p)/(2-7p)2. Equation 
(jl|) then implies 

4 

m* = (/'*ecorc — j=kpA^''^ — , (2) 
vrv-J ^ 

which is a generalization of the isothermal accretion rate to the nonthermal case. 

For the polytropic sphere, p oc r"^'' and M oc r^~^p^ which implies p oc 
^-fcp/(3-fcp) ^ ^-2/(4-37p)_ Thus we have oc m^p^/^ oc mi"^'^^^"^^''^ for the 
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mass dependence of the accretion rate. Integration yields 

4-37P 

(3) 

where is the final stellar mass, which is attained at a time t*/ (McLaughlin 
& Pudritz 1997, hereafter MP97). Note that for 7^ < 1 the accretion rate 
accelerates (MP97). Equation (^) implies 




m* . . „ X m*/ / t 



m. = (4-37p)-l = (4-37p)-^ — . (4) 

t t^f \t^f J 

As discussed by MP97, termination of the accretion breaks the self-similarity 
once the expansion wave reaches m^j. This occurs at a time they label te^f,, 
which is about 0.45t*j. Thereafter, equation (^) becomes approximate, but 
they argue that the approximation should be reasonably good. From equations 
(|l|) and @), the star-formation time is 

t.f = (5) 
Using the results of MP97, we can evaluate (p^, in the non-magnetic case, 

■(2-7p)2(4-37p)(^-'^>)/2mo^'/^'-'^^^ 



(4-37p)^=vrV3 



8(5-37p)/2 



(6) 



where tuq is tabulated by MP97. For example, for the singular isothermal sphere 
the parameter mo = 0.975 and 4>*non = 0.9757r\/3/8 = 0.663. For other values 
of 7p in the range < 7^ < 1, </>*non — 1.13/(1 + 0.77^). We conclude that 

— m^:/tff to within a factor 1.5 for spherical cores in which the effective 
sound speed increases outward. 

We can estimate the effect of magnetic fields on the accretion rate from 
the work of Li & Shu (1997), who considered collapse of self-similar, isother- 
mal, magnetized, toroidal clouds. The equilibrium surface density is S = 
(1 + Hq)c^^/{ttGw), where w is the cylindrical radius and Hq is a parame- 
ter that increases from zero as the magnetic field is increased. They show 
that the accretion rate is = 1.05(1 + Hq)c\-^/G, which is larger than the 
isothermal case by about a factor (1 + Hq). However, equation (||) predicts 
oc 0^ ^3/2/^3/2 ^ 5^3/2^3/2 ^ (I Hofl"^ . To reconcile this result 

with the correct value, we require 0* ~ i?5>*non/(l + Hq)^/"^. They regard i^o ~ 1 
as typical for low-mass star formation, which would lead to a reduction in the 
accretion rate (as expressed in eqs. [|| and [^) by a factor of about 1.4. 

The power-law structure of cloud cores is truncated by the ambient pressure 
in the star-forming clump, which is thus equivalent to the core surface pressure, 
Ps- We parameterize the rate of individual core collapse in terms of this pressure, 
which we then relate to the clump's surface density. A higher ambient pressure 
means that a core of a given mass is truncated at higher density and thus its 
free-fall and star formation timescales are shorter. 
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Since ps = Ps/cl, M = {kpc^r)/G and p = A(? j {lixGr^) , the density at 
the surface of the core is 



1/4 

(7) 



The general expression for the accretion rate in terms of the surface pressure 
and the final stellar mass can now be inferred from equations (|^) and (^): 



3(2-fcp) 

^ ' - i / \ 2{3-fcp) 



(8) 



10/ \ iU"n.cm " / \ rn^f 
and the corresponding value of the star-formation time is 

3/8 



^4 4 - 37p \ / m*/ ( 10^ K cm-^' 



The accretion rate and star formation time depend to some extent on the 
value of kp in the protostellar cores. No data are available on the structure of 
cores that are forming very massive stars. For "high-mass" cores in Orion, Caselli 
& Myers (1995) find kp ~ 1.45 with a dispersion of ±0.2. According to van der 
Tak et al. (2000), the clumps in which high-mass cores are embedded have 
values of kp ranging from 1 to 2, centered around 1.5. For numerical evaluation 
of the accretion rate, we shall adopt kp = 1.5, which corresponds to 7p = 2/3, 
kp = 1, A = 3/4, and 4> 

*non — 0.90. In this case, the protostellar accretion rate 
varies linearly with time, oc t. So long as ^ oc 7p = 2{kp — l)/kp is not 
near zero, the star-formation time is relatively insensitive to kp. For kp = 1.5, 
the numerical coefficient in equation (^ is 3.49 x 10~^/(1 + Hq)^/"^ Mq yr~^ and 
m=K oc {m^/m^,f)^'^ . Similarly, the numerical coefficient in equation (|9|) becomes 

1.72 X 10^(1 + i?o)^/^ yr. 

At present, direct observations of the total pressure at the surface of a core 
that will form a massive star are not available, so we must estimate this pressure 
theoretically. Since the core is embedded in a turbulent clump, its pressure will 
fluctuate, but the time-averaged surface pressure for the typical core should be 
approximately the same as the spatial average of the pressure in the clump, 
Ps = <ppP, where <j)p is a. parameter of order unity. The mean pressure in a 
clump is (MT02): 



3vr/gas</'Bavir 



20 



GS^ (10) 



tvir 



where /gas is the fraction of the clump mass in gas (as opposed to stars), 
5{a'^)Rci/GMci = Mvir/Mci is the virial parameter, a is the one-dimensional 
velocity dispersion, (x) indicates a mass-averaged quantity, (j)B = (c^) / {cr'^) — 
1.3 + 1.5/m^ represents the contribution of the magnetic pressure, mA is an 
average of the Alfven Mach number, and S = Mci/{itR'^i) is the surface density 
of the clump, including the mass of any embedded stars. In terms of Svir, the 
surface density determined from the virial theorem, we have E = Svir/ovir- 
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^; the corresponding 
r X 214) mag. These 



Regions of high-mass star formation studied by Plume et al. (1997, hereafter 
P97) are characterized by virial masses Mvir ~ 3800 Mq and radii ~ 0.5 pc. 
The corresponding mean column density is Y^^ir — 1 g cm~ 
visual extinction is Ay = {Nu/2 x 10^^ cm~^) mag = (Svj 
column densities are far greater than those of giant molecular clouds (GMCs) 
(0.035 g cm~^ — Solomon et al. 1987) or of regions of low-mass star formation 
(the average S in the C^^O cores in Taurus is 0.032 g cm~^ — Onishi et al. 1996). 
On the other hand, the central stellar surface density in the Orion Nebula Cluster 
is about 1 g cm~'^ (Hillenbrand & Hartmann 1999). The mean density in the 
Plume et al. clumps is tt-h — 2 x 10^ cm^^. As discussed by MT02, the virial 
parameter Ovir — 1-3 — 1.4 for GMCs, whereas it is quite close to unity for cores 
that are actively forming stars. Although there is no direct determination of 
avir for massive star-forming clumps, P97 regard the virial mass as the most 
accurate, and we adopt Ovir = 1 for our numerical estimates. 

MT02 discuss the estimates for the other parameters, (pB — 2.8, /gas — 2/3 
and (^p ~ 2 needed to determine P,. Thus from equation (10) 



Ps = 2.28 X 10' 



' (/gas^P'^Bttvir) 5]^ K Cm" 



8.5 X 10^ 



K cm 



-3 



(11] 



With ec 



0.5 and Hq = 1, the corresponding protostellar accretion rate is 



4.75 X 10-^e,^ 



/4 {fgas<PpaviT(pBf^^ f m^f 



4.6 X 10" 



(l + Fo)i/2 V3OM0 

0.5 



3/4 



5.3/4 



0.5 



m 



30 Mq 



Mq yr- 



Moyr-^ 
(12) 



and the star formation time is 
1.26 X 10^ 



tcore 



Uga! 

1.29 X 10^ 



m 



*f 



30 Mp 







1/4 



^3/8 
(t>B 

5.-3/4 



m 



*f 



30 M, 



o 



1/4 



yr. 



S^3/4 (13) 
(14) 



For example, with the fiducial values of our parameters, a 100 Mq star forming 



in a clump with S = 1 g cm^ has a final accretion rate of 1.1 x 10^ Mq yr~ 
and a star-formation time of 1.75 x 10^ yr. Our result for the accretion rate is 
somewhat lower than that of Osorio et al. (1999), who found accretion rates 
by matching the infrared spectra of embedded protostars. Our result for the 
star formation time is comparable to that of Behrend &: Maeder (2001), who 
performed an analysis based on protostellar outflows, although our result for 
the time evolution of the accretion rate is quite different. 

The above solutions apply to purely nonthermal cores. We have also al- 
lowed for an additional thermal (T ~ 100 K) component of pressure support 
(MT02). For typical high pressure star- forming regions cores with M ^ 1 Mq 
are dominated by the nonthermal component. This implies that cores are super- 
sonically turbulent and therefore clumpy, which means we expect the accretion 
rate to exhibit large fluctuations about the mean value. Radiation pressure can 
reduce the accretion flow to the most massive stars from the values predicted 
above (Wolfire & Cassinelh 1987; Jijina & Adams 1996; Tan & McKee 2002b). 
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-10 12 

Figure 1. Evolution of protostellar and L* in a S = 1 g cm~^ 
clump. For references on HMC luminosities see MT02. 

3. Protostellar Evolution, Luminosities and Outflows 

The properties of accreting protostars depend on their accretion rates (e.g., 
Stahler, Shu & Taam 1980; Stahler 1988; Palla & Stahler 1992; Nakano et al. 
2000). In particular, the accretion luminosity Lace = /acc Cm* m^</r*, where /acc 
is a factor of order unity accounting for energy radiated by an accretion disk 
or used to drive protostellar outflows, and the stellar radius may depend 
sensitively on m^:. Massive stars join the main sequence during their accretion 
phase at a mass that also depends on the accretion rate (Palla & Stahler 1992). 
The intensity of protostellar outflows, at least from low-mass stars, is thought 
to depend on the accretion rate and the Keplerian velocity, vk, at the stellar 
surface (e.g. Shu et al. 1994). To predict protostellar properties and thus 
compare theory with observation, we require a model of protostellar evolution. 
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As previous investigations have mostly considered constant accretion rates, 
we have developed a simple model for protostellar evolution based on that of 
Nakano et al. (2000), which allows easy implementation of time dependent 
accretion. The model accounts for the total energy of the protostar as it accretes 
and dissociates matter, expands and contracts, and, if the central temperature 
Tc ^ 10^ K, burns deuterium. We have modified Nakano et al.'s model to 
include additional processes, such as deuterium shell burning, and calibrated 
these modifications against the more detailed calculations of Stahler (1988) and 
Palla & Stahler (1992). Protostellar evolution is followed from m* = 0.1 Mq to 
the end of accretion (m* = m*j). Massive protostars reach the main sequence 
(Schaller et al. 1992) before they have finished accreting. The evolution of 
protostellar radius as a star grows in mass is shown in Figure 1 for several 
different final stellar masses, all forming in a clump with S = 1 g cm"^ and 
T = 100 K. Also shown is the evolution in bolometric luminosities (including 
the contribution from an accretion disk) , and our use of these models to predict 
the masses and accretion rates of protostars thought to be illuminating several 
nearby hot molecular cores (HMCs). MT02 consider the effect of varying kp and 
S on these estimates and also compare results with Osorio et al. (1999), who 
have modeled the same sources. 

Powerful outflows are ubiquitous from regions of both low and high-mass 
star formation (e.g., Richer et al. 2000). Models of magnetically driven outflows 
from a disk (Blandford & Payne 1982) and from the inner edge of a disk (Shu 
et al. 1994) have been presented. A common feature of these models is the 
production of a bipolar outflow with momentum distribution oc (sin^)"^ for 
9 > Oq, where 9 is measured from the outflow axis and 9o ~ 10"^ (Matzner & 
McKee 1999). On scales large compared to the source (Matzner & McKee 1999) 

dpw Pw 

~ 47rln(2/^o)(l + ^0 -cos2^)' ^ ^ 

We parameterize the total momentum flux escaping in the outflow via 

Pw — fw,esc^wVw — fw,escfw^*Vw — 4'w'^*VK i (1^) 

where fw,esc is the fraction of outflow momentum that escapes from the core, 
= Tfiyj/m^, = fw,escfwVw/vK and vk is the Keplerian velocity at the 
equatorial radius of the star. Najita & Shu (1994) considered the acceleration of 
winds from an accreting 0.5 M© protostar. For several different X-wind model 
boundary conditions their results gave quite a large variation in ~ 0.1 — 0.8, 
but only a small spread in fwVw/vx — 0.6 within ~ 30%. Najita & Shu (1994) 
also showed that Vyj is approximately independent of 9 so that rriyj has the same 
distribution as p^, and thus /w,esc applies to wind mass as well as momentum. 
If all the protostellar wind material escapes in two cones of solid angle ^w,esc 
(opening angle 9w,esc), and if core material originally inside these cones is also 
swept up, then 

_ 1 - 20^,esc/(47r) 

Ccore ~ 1 I f f • y ' I 

-L ~r Jw,escjw 

For ecore = 0.5 (Matzner & McKee 2000) and fw = 0.1,0.3,0.8 we solve equa- 
tions (ll|) (with rhyj) and (0) for f^^^sc and 9^^esc- We find /^,esc = 
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0.89,0.86,0.79 and 

Gw,esc — 57°, 51°, 36°, respectively. The escape efficiencies 
are quite higli because of the concentration of the wind near the outflow axis 
(eq. [|l5|). Assuming the wind is atomic, the swept-up gas is molecular and they 
escape together in a momentum conserving shell, the fraction of momentum 
leaving the core that is atomic in these three cases is 0.49, 0.62, 0.78, respec- 
tively. 

We adopt = 0.3 and /«,,esc = 0.86 as our fiducial values so that cp^ ~ 0.5, 
since fwVw/vK — 0.6 (Najita & Shu 1994). This estimate is broadly consistent 
with observations of outflows from low-mass stars (Bontemps et al. 1996). We 
apply the above model for low-mass protostellar outflows to stars of all masses 
as they accrete. This application is uncertain: for the X-wind model (Shu et al. 
1988) to apply, the stellar surface field strength required by a 10 M.q protostar 
with radius lOR© accreting at a few xlO~^M0 yr~^ is ~ 10'^ G. The momentum 
injection rates and total momenta from protostars forming in a S = 1 g cm~^ 
environment are shown in Figure 2. A fit to these results is 

P«,-Pi<J=68(^) (^^'^Mekms-i. (18) 

This relation depends only weakly on S (TM02). 



4. Probing the Rate of Star Formation in Clusters with Outflovi^s 

We model the formation of star clusters with M^:f = M^f /m^,f stars, where M^f 
is the total stellar mass and m^,j is the mean individual stellar mass. Stars are 
drawn from a power-law initial mass function (IMF) truncated above and 
below m^i and of the form dF/d In m^, = A^,m~°^* , where F is the fraction of all 
stars that have a mass greater than m*. The Salpeter (1955) IMF has a* = 1.35, 
= 0.0603 and = 0.353 M© for m^i = 0.1 M© and = 120 Mq. We 
adopt this as our fiducial IMF. While observational determinations of the IMF 
indicate deviations from this simple power-law form at lower masses, these have 
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little effect on most of the feedback processes operating in a protocluster, which 
are controlled by the massive stars. The total protostellar wind momentum, 
Pw,tot) from the forming star cluster does receive important, approximately equal, 
contributions from all decades of the IMF, since a* — ctp. However, even in this 
case, the variation caused by uncertainties in the low-mass IMF is minor. The 
low-mass stars do dominate the total stellar mass, and this may be parameterized 
via /i/i, the total stellar mass per high-mass (m* > 8 Mq) star, most of which 
end their life in a core-collapse supernova explosion. For our adopted IMF, 
^j,h = 134 Mq, which is close to the value for the Miller & Scalo (1979) IMF. 

We use Monte Carlo simulations to find the median values of the post- 
accretion luminosity, L^j, and total protostellar outflow momentum released 
during formation of a cluster of mass M^j. Typical lower- mass clusters fail to 
fully sample the IMF. For M^j/Mq ^ 1600, L^f/M^j ~ 630 L© Mg^ while for 
100 ^ M^j/Mq ^ 1600, L^f/M^f ~ 436(M*//lOOOM0)O-8 Lq Mq^ For the 
outflows, with the flducial values of the parameters S = 1 g cm^^, T = 100 K, 
kp = 1.5 and (p^ = 0.5, we find 



For M^f > 1000 Mq, Pw,tot/M^f ~ 87(0^/0.5) kms"^. The loss of momentum 
in outfiow-outfiow interactions would lower these estimates. However, because 
of the collimation of momentum in jets (eq.||l^), we expect this effect to be 
relatively minor, particularly for smaller clusters with fewer simultaneous inde- 
pendent outfiows. 

Unless triggered by an external infiuence such as a supernova blast wave, 
the overall rate at which stars form within a clump is not expected to exceed 
~ M*j/tg, where % = {3tt/32G'P)^^'^ is the free-fall time of the clump evaluated 
at the mean density (stars + gas) at a typical stage in the evolution (i.e. when 
half the stars have formed). We write the star formation rate in terms of the 
numerical parameter rj > 1 as 

r]ts rits 3rjtf[ 

where M^^^typ = 0.5M^f is the typical current mass of stars, Md^typ is the typical 
current total mass of the clump, and, for ecore = 0.5, M^^^typ = -^ci,typ/3, if 
the protostellar outfiows are efficiently ejected from the clump and most of the 
clump mass is in cores (see TM02). The star formation rate is then 

M. = 9.7 X 10-=^ ( -] ( -^^^y^^ Mq yr-i (21) 



\r]J VIOOOMq 
and the cluster formation time is 



^cluster = f?iff = 6.86 X 10^77 ( ^"^'^/^ ^ ^ S"^/^ yr. (22) 
' ' V 1000 MqJ ^ ^ ' 



We wish to estimate the cluster star formation rate, and thus t/, from ob- 
servations. One method is to relate the summed total of many individual stellar 
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accretion events to the production of a protocluster wind, which is the super- 



position of many individual bipolar protostehar outflows. From equations (19) 



and (^2]) the expected momentum injection rate from a forming star cluster is 



^cluster 

~\oI,)\'Wm^) ^ MoJtms yr < 15UU 



s/4 (23) 
f (fe) (wit) ^^/^ M« km s- yr-i > 1500. 

We compare equation (^) against observations of outflows from protoclus- 
ters (Choi et al. 1993, hereafter C93). The observational determination of p^^j^tot 
involves measuring the total outflow momentum over some region and identify- 
ing a characteristic timescale of the flow over the scale of this region. The CO 
observations of C93 cover regions with sizes similar to Rc\ for which virial mass 
estimates have been made from CS transitions, tracing dense gas. The spatial 
extent of the outflow observations is limited by sensitivity and the flows almost 
certainly continue beyond Rd. This is also evident from the fact that the flow 
timescales are much shorter than These considerations allow us to neglect 
the effects of evolution in determining an estimate of the instantaneous star for- 
mation rate from the observed flows. Though poorly resolved, the morphologies 
of the flows appear quasi-spherical, as would be expected if there were multiple 
driving sources. 

From CO line strengths, C93 derive the mass of the flow, correcting for opti- 
cal depth and excitation temperature effects and assuming standard abundance 
ratios. The line of sight momentum is then summed over velocity channels. We 
have further corrected these estimates for projection effects (a factor of two in- 
crease for a spherical flow). The total momentum (and mass) of the flow may 
be greater by the amount contained in atomic gas. For example we have seen 
in our fiducial model for protostellar winds from an individual core that 62% 
of the total momentum is in this form (i.e. resides in gas that was processed 
through the accretion disk). As this flow sweeps up gas beyond the core, the 
atomic fraction is reduced. TM02 estimate that about 25% of the total momen- 
tum is in atomic form, so to recover the total momentum we therefore increase 
the observed value by 33%. This question may be examined observationally, 
though current surveys have placed only upper limits on the fraction of atomic 
momentum in the flows (Choi et al. 1994). 

Given estimates of mass and momentum, we determine a characteristic flow 
velocity, Vp, and thus an outflow timescale tp = Rp/vp (here we use the data 
uncorrected for projection effects, which should approximately cancel in this 
ratio). The value of pu;,tot is then calculated, and with knowledge of the clump's 
mass and surface density, used to estimate r/. The properties of our sample of 
protoclusters and their outflows are listed in Tables |l] and |2[ 

Excluding the cases of Cep A and NGC6334 I(N) (below), we find the loga- 
rithmic mean value of r/ is 10.4. The standard deviation is about a factor of two, 
consistent with the observational uncertainties. This gives a cluster formation 
time (eq.Q) of 7.1 x lO^(Mci,typ/lOOOM0)i/^S-3/'^yr, which is consistent with 
estimates of ^ 10^ yr from the spread in pre-main-sequence ages in the Orion 
Nebula Cluster (Palla k Stabler 1999), which has Md ~ 4500 (Hillenbrand 
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Table 1 . Properties of Protoclusters"^ 



Name 


d 




R 


S 






(kpc) 


(Mo) 


(pc) 


(g cm^^) 


(105 Lq) 


S140(IRS1) 


0.9 


128 


0.14 


0.43 


0.05 


GL490 


0.9 


91-151 


0.06-0.083 


1.68-1.44 


0.017-0.024 


GL2591 


1.0 


268-320 


0.13-0.11 


1.05-1.76 


0.225 


NGC2071 


0.39 


484 


0.69 


0.068 


0.0052 


W3(Main)'' 


2.3 


610 


0.35 


0.33 


5.2 


Cep A 


0.73 


787 


0.14 


2.67 


0.17 


W28(A2) 


2.0 


982 


0.18 


2.0 


1.3-1.9 


NGC6334 I(N) 


1.7 


~1700 


0.37 


0.83 


< 0.1 



Data and lists of primary references for all sources except NGC6334 I(N) (Megeath & 
Tieftrunk 1999) are in C93, P97, and van der Tak et al. (2000). 

W3(Main) appears to be a relatively evolved system, containing numerous H II regions, so 
the virial mass estimates based on observations (P97) of the gas today probably underestimate 
the initial mass. This can explain why Mci/3 is an insufficient stellar mass to account for Lboi- 



& Hartmann 1998). We note that the effect of wind-wind interactions, which 
we have ignored in our analysis, would shorten this estimate by the fraction of 
wind momentum that is dissipated. 

Cep A is the most extreme of the C93 protoclusters in terms of its surface 
density. The estimates of virial mass and surface density are made uncertain by 
the fact that 2D maps of the dense gas show moderate elongation (Y. Shirley, 
private communication). Furthermore, the outflow observations do not symmet- 
rically probe the central region of Cep A (M. Choi, private communication) and 
so the total Pw,tot may be underestimated, thus artificially raising r/. For these 
reasons we have excluded Cep A from our statistical sample. NGC6334 I(N)'s 
outflow was observed by Megeath &: Tieftrunk (1999) and the outflow parameters 
determined using different techniques to those of C93. This source is thought 
to be in an extremely early stage of evolution; for example the observed outflow 
extends over only a small fraction of the clump, which has a relatively low lu- 
minosity. Again we exclude this source from the statistical sample, but consider 
it in the context of an example of the earliest stage of cluster formation. 

A common observational diagram used in the study of outflows from proto- 
stars and protoclusters plots Pw,tot versus Lboi- We consider a simple description 
of the evolution of protoclusters in this diagram. Assuming a constant star for- 
mation rate (eq.|21]) implies that, averaged over a population of clusters, Pi„,tot 



is a constant over most of the evolution. In the earliest stages the luminosity is 
small and dominated by that resulting from accretion. The minimum luminosity 
is ~ Gm=KM*/r*(m=K), where ^ m^f is the mean protostellar mass. We have 

-M* = Pw,tot/ {4'wVk), where vk is the Keplerian velocity at the stellar surface 
of an average protostar. The mean radius is about 3 R© (Figure 1) and thus 

L„,. = 4.3 X 10' (f) (^^) L _. ) Le. (24) 

V0«,/ \0.3Mq r^{m^) J VMokms ^yr V 
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Table 2. 



Tan & McKee 
Observed Outflows from Protoclusters" 



Name 


Rp 


tp 


Pw, tot 


Pu),tot 






HV(EHV)'' 


HV(EHV) 


HV(EHV) 


HV(EHV)[TOT]'= 


V 




(pc) 


(103 yr) 


(M0kms-l) 


(lO-3M0kms-iyr"^) 




S140(IRS1) 


0.20 (0.15) 


31 (6.6) 


307 (10.2) 


10.1 (1.6) [15.4] 


4.4 


GL490 


0.087 (0.087) 


9.7 (2.4) 


60 (11.8) 


6.2 (4.9) [14.8] 


9.5-13.2 


GL2591 


0.097 (0.055) 


11.9 (2.4) 


131 (3.11) 


11.0 (1.3) [16.4] 


15.7-27 


NGC2071 


0.051 (0.044) 


4.2 (1.5) 


13.2 (3.54) 


3.1 (2.4) [7.3] 


7.7 


W3(Main) 


0.39 (0.20) 


29 (7.5) 


204 (16.2) 


7.1 (2.2) [12.2] 


18.4 


Cep A 


0.11 (0.11) 


11.3 (3.0) 


156 (15.3) 


13.8 (5.1) [25.1] 


54 


W28(A2) 


0.14 (0.14) 


10 (3.4) 


889 (82.7) 


88 (24) [149] 


8.9 


NGC6334 I(N) 


0.22 (n/a) 


9.5 (n/a) 


312 (n/a) 


33 (n/a) [44] 


25 



" Data from C93, except for NGC6334 I(N) (Megeath & Tieftrunk 1999). 
C93 outflow data are divided into high (HV) and extremely high (EHV) velocity parts. 
Includes correction factor of 1.33 for atomic gas (TM02). 



Observed protoclusters, selected for showing signs of massive star formation, 
will tend to be over- luminous compared to equation (p4|). Those in the earliest 
stages of evolution should lie systematically closer to the birthline. 

The maximum luminosity is approximately the post-accretion luminosity, 
L^f, for all but the smallest clusters (where accretion can make a significant 
contribution). We have determined L=kj(M*j) from Monte Carlo simulations 
(above). Relating M*j = 2Mci^typ/3 for a typical forming cluster (TIVI02), we 
then eliminate Md^typ in equation (p^). The resulting expression defines the 
"cluster birthline" (the 3 cases correspond to values of M^j/Mq of 10^ — 10^, 
10^ — 1.6 X 10'^ and > 1.6 x 10'^, respectively): 



L 



*f 



3.5 X 1070^2.05.-1.5 (^_^^^^^-^ ■ X 103 ^ ^ ^ 4.4 x 10 

7.9 X lo70^2-^S-i-« ( M0kts-VO '''^Q X ^ ^ 1-0 X 10' 



1.1 X 1O70-/^S- ( mJs^V.-0 1.0 X 10^ ^ 

(25) 

where (Pl = (</>^/0.5)(10/r/). 

The protocluster and cluster birthlines are shown in Figure 3. In Figure 3a 
we also plot the results of Monte Carlo simulations of many clusters for two 
cases: = 200,1000 Mq. For each cluster, stars are drawn randomly from 

the IMF and stochastically in time and then followed through our protostellar 
evolution model for pressures corresponding to S = Ig cm~-^. The central/outer 
lines show the median/68-percentile values of each cluster distribution. Points 
mark equal time intervals in the evolution. In Figure 3b we plot the observed 
protoclusters, giving pw,tot a factor of four uncertainty. Sources in an early 
stage of formation (NGC2071 and NGC6334 I(N)) are close to the protocluster 
birthline, while those in a late stage (W3(H2 0) — as evidenced by its multiple 
H II regions) are close to or below the cluster birthline. 
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234567 234567 



Figure 3. p«,,tot versus Lboi- (a) Left: Theoretical birthlines and 
results of Monte Carlo simulations, (b) Right: Observed protoclusters. 



5. Conclusions 

We have presented a theoretical model for high-mass star formation that deter- 
mines the accretion rate in terms of the final and instantaneous masses of the 
star, the ambient pressure surrounding the gas core, and the polytropic index of 
pressure support. We related the pressure due to self-gravity in the centers of 
typical Galactic star- forming clumps to their surface density; P/k approaches 
10^ K cm~^ for typical values of S ~ 1 g cm~^. These pressures allow a 100 Mq 
star to form in ~ 10^ yr with a final accretion rate ^ 10^'^ Mq yr^-*^. Model- 
ing protostellar evolution, we predicted the properties of several nearby massive 
protostars. Under the hypothesis that the outflows from massive protostars are 
driven magneto-centrifugally in the same manner as those from low-mass proto- 
stars, we found that each decade of the IMF contributes approximately equally 
to the total outflow momentum. This helps to explain the generally poor overall 
collimation of outflows from high- mass star- forming regions. We utilized these 
models and observations of outflow intensities from several nearby protoclusters 
to estimate their star formation rates, finding that these clusters are typically 
forming over at least several free-fall times. We have proposed theoretical birth- 
lines in the Pw,tot versus Lboi diagram to delimit the phase of cluster formation. 
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